Construct order parameter from the spectra of mutual information 
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In this paper, we try to establish a connection between a quantum information concept, i.e. the 
mutual information, and the conventional order parameter in condensed matter physics. We show 
that a non- vanishing mutual information at a long distance means the existence of long-range order. 
By analyzing the entanglement spectra of the reduced density matrix that are used to calculate the 
mutual information, we show how to hnd the local order operator used to identify various phases 
with long-rang order. 
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I. INTRODUCTION 



At absolute zero temperature, a quantum many-body 
system may undergo a phase transition as the system's 
parameter varies across a critical point Except for 
some unconventional phase transitions, continuous quan- 
tum phase transitions can be well characterized in the 
framework of the Landau's symmetry-breaking theory. 
In the Landau's theory, the concept of the order pa- 
rameter plays a central role. To characterize a given 
ordered phase, people usually introduce an order parame- 
ter which is nonzero in the symmetry-breaking phase, but 
vanishes in the other phases. In the recent years, much 
attention has been paid to using quantum information 
approaches to investigate quantum phase transitions. 
A typical example is the role of entanglement in quan- 
tum critical phenomena Since the quantum phase 
transitions occur at the absolute zero temperature, they 
are purely driven by quantum fluctuations. People be- 
lieve that the quantum entanglement, as a kind of quan- 
tum correlation, should play a very important role in the 
quantum phase transitions. This point has been proven 
to be true in the later hundreds of works on the issue(for 
a review, see [9| and the references therein). The entan- 
glement manifests interesting properties, such as scaling 
[3, |4j, singularity, or maximum [5|,|6| etc, in various quan- 
tum phase transitions. 

However, it seems to us that, besides its interesting 
behavior around the critical point, the entanglement ac- 
tually cannot tell us much information about the corre- 
sponding phase itself. That is, while the entanglement 
helps us to witness the quantum phase transition, a non- 
trivial question is how can we learn the corresponding 
order parameter of the ordered phases from the behav- 
ior of the entanglement. This question leads to the main 
motivation of the present work. We will show that the 
non-vanishing mutual information Q at a long distance 
not only means the existence of long-range order, but also 
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can tell us the potential order parameter if we analyze the 
entanglement spectra of the corresponding reduced den- 
sity matrices. A possible scheme to construct the order 
parameter will be discussed in detail. 

On the other hand, the order parameter plays a very 
important role in condensed matter physics. In order to 
study the ground-state properties of a quantum many- 
body system, people usually assume a special order pa- 
rameter in various methods. For instance, in the dynamic 
mean-field theory, a suitable local order parameter usu- 
ally is assumed at the very beginning. However, such 
an order parameter strongly depends on the physical in- 
tuition of the investigator. That means, if the guessed 
order parameter is wrong, then all analysis based on the 
assumption of the order parameter are not physical or 
even wrong. Therefore, a scheme that can tell us how 
to find an appropriate order parameter is instructive for 
many condensed matter theorists in their studies. 

The paper is organized as follows. In Sec. II, we firstly 
discuss the relation between the existence of long-range 
order and the non- vanishing behavior of the mutual in- 
formation. In Sec. Ill, we propose a possible scheme to 
find the local order operators by analyzing the entangle- 
ment spectra of the mutual information. Three simple 
applications are given in the same section. Finally, we 
give a discussion and a brief summary in section IV. 



II. NON- VANISHING MUTUAL 
INFORMATION AND LONG-RANGE ORDER 



Let us consider a general quantum many-body system 
described by the Hamiltonian H = hi, where hi is the 
local Hamiltonian. hi can be the kinetic energy, interac- 
tion, etc. To learn the potential long-range correlations 
existing in the system's ground state we focus on 
two small blocks, i and j, separated by a distance \i — j\. 
Without loss of generality, we assume that the block has 
a minimum size with which the block can capture the 
potential long-range correlation existing in the system, 
and the distance \i — j\ is comparable to the system's 
size. The reduced density matrix of the block i can be 
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expressed as 



in the basis of 



'Hi/ 1 



Insert the identity 



(p'\Pi\p) = tr(oi/i'PoL)> 



(1) 



where p = \*5>) (\&|, tr is used to trace out all other de- 
grees of freedom of the system except for the block i, ai^ 
is the annihilation operator for the state \p) localized at 
the block i. The operators a.^s satisfy the commutation 
(or anti-commutation) relation for bosonic (or fermionic) 
states. Similarly we can define the reduced density ma- 
trix of the two blocks i and j as 

{p'v'\p iuj \p,v) =tx(a ill 'a j u'pa] v a\ tl ). (2) 

The reduced density matrices are positive semidefinite 
and can be normalized as 

ti(pi) = ti(pj) = tr(piuj) = 1. (3) 

So if we diagonalize these matrices 

]p^\^)(^\, (4) 



V 

PiUj — ^ Qfj,u\<Piii,jv) ((f>in,ju\ 



(5) 
(6) 



for a given reduced density matrix, say pi, the diagonal 
elements {p^} then have a probability interpretation in 
the corresponding eigenstate space. The von-Neumann 
entropy 



A* 



Pfi log 2 Pp. 



(7) 



measures the entanglement between the block i and the 
rest of the system. 

To study the correlation between the block i and j, we 
introduce the mutual information 



= S{pi) + S (pj) - S (piuj) 



(8) 



In the quantum information science, the mutual infor- 
mation measures the total correlation between the two 
blocks. It has also been used to study quantum critical 
phenomena in recent years [7, 8]. For the completeness 
of the work, we show that the non- vanishing behavior of 
the mutual information at a long distance leads to the 
existence of long-range order, as discussed in Ref. Q. 
For this purpose, we express the mutual information as 

S(i\j) = tr(p iUj log 2 pujj) - tr(pju f log 2 Pi ® Pj)- (9) 

Then 

-^(fawvlpiuj log 2 (p* ® Pj)\^,ju), (10) 



Yluv IwwXwyVil — 1) the mutual information 
becomes 

s (i\j) = £ q » v log ' 2 q » v ~ P »'"',n>>qn'i>' log 2 (p M P^) 

(11) 

with 



P, 



<pjv') (<PW tPjv>\<f>i»,ju) (12) 



which satisfies 

P)x'u'^[lV ^ 0; ^ P fl' v' ,flU — lj ^ PjX 1 U 1 .[lV 1- (13) 

H'v' pv 

Eq. (ITT1) can be rewritten as 

= E q » v lo §2 - E ^W'"' l °g2(Pn'Pv') 

(14) 

Now we divide the problem into the following two cases. 

Case I): For any matrix element in pujj, if we can al- 
ways find (a\^ai^a] u aj vl ) = (at t Oi^/)(ot l/ o 3V ) such that 
Piuj — Pi ® Pj, then P^'u'^u is a unit matrix and we 
have S(i\j) — 0. In this case, any linear superposition 
of a| aj^/ is not correlated, so the system does not have 
any long-range order. 

Case II): If there exist some elements that 
(<4n a in' a jv a jv) ^ (4^')( a ju a jw), then p iuj ^ pi®pj, 
Pfj,'v',fa/ in Eq. (fT2")l is no longer a unit matrix. Since the 
log function is concave, 

^P^v',pv^Og 2 {P^Pv') < l0g 2 ^P^v'.pvPpPv 

we have 

S{i\j) > 2l<liJ,v\og 2 Q lxv 

fa/ 

wherewhere 



So the mutual information is non-vanishing 
S(i\j) > 0. 



(16) 



In this case, pujj 7^ Pi ® Pj- Actually, S(i\j) > is the 
sufficient and necessary condition for the existence of the 
correlation between two blocks. That is, if and only if 
S(i\j) > 0, we can find matrix elements satisfying 

(a\^a itl >a] v a jvl ) ^ (a^a^at^a^). (17) 



3 



Then (aj a^') plays a role of local order operator at site 
i. Nevertheless, there are usually many terms satisfying 
the above inequality, the operator (a\ di^) is not a good 
candidate for long-range order parameter. 

To find a good candidate of the order parameter, we 
need to analyze spectra of the mutual information. To 
proceed, we would like to mention some properties of the 
mutual information: 

i) If the rank of pi is 1, then there is no correlation 
between pi and pj. 

Prove: If the rank of p, is 1, then pi is a pure state 
and S(pi) = 0. Since S(i\j) = S (pi) + S (pj) - S (p iuj ) 
is non-negative, so S(i\j) = 0. There is no correlation 
between the blocks i and j. 

ii) If the rank of pi is £, for any p > £, the operator 
a \^ a i^ docs not correlated. 

Prove: If 

Pi\<Pip)=0, (18) 

Then 

- 0. (19) 
The correlation function becomes 

= ( fl L a v I W Pj<" ) (<*V Pj"' I a ]n a m) 

\x'v' 

= 0. (20) 

iii) For p ^ v, (a^a^) = (a\ v a itI ) = 0. If 

a l A i a ^ a ] A1 a jv) 7^ 0' a ln a ™ + a L a »M is a g 00 ^ candidate 
of the order parameter. 

These three properties will be used in the next section 
to construct local order operator. 

III. CONSTRUCT LOCAL ORDER OPERATOR 

In this section, we show how to construct the potential 
local order parameter from the spectra of the reduced 
density matrices. We have show that the long-range or- 
der exists if and only if the mutual information does not 
vanish at long distance. However, the mutual informa- 
tion can not tell us directly which kind of order exists 
in the phase. To find the potential order parameter, we 
need to look into the spectra of the mutual information. 

At the beginning, we need to decide the minimum size 
of the two correlated blocks. If the block size is too large, 
even if the mutual information does not vanish, the final 
order parameter obtained might be too complicated to be 
well understood. On the other hand, if the block size is 
too small, the corresponding mutual information might 
not be able to find the order parameter. For instance, in 
the spin dimer or trimer ordered phase, the mutual infor- 
mation between two single spins usually decays exponen- 
tially because the spin-spin correlation function decays 



exponentially too. Therefore, to decide an appropriate 
block size is important. For this purpose, we can start 
from a relatively large block size. In this case, it is easier 
to witness if there exist long-range order in the system. 
Then we reduce the size of the blocks step by step until 
find the minimum size of the correlated blocks. Once the 
minimum block size is decided, we continue to find the 
potential order parameters from the mutual information 
of the two blocks. 

Diagonal long-range order: Usually, there are two 
kinds of long-range order. One is the diagonal long-range 
order, the another is the off-diagonal long-range order. 
Whether the order is diagonal or off-diagonal depends 
on the basis used to define the reduced density matrices. 
Without loss of generality, we use the eigenstates |w) 
of pi to define the local modes. Then the set {pi^} is the 
probability distribution of these modes. 

We now consider the first case of diagonal long-range 
order. In this case, if we express the combined block p iU j 
in the basis l^Pi^fju), we can find that there exist some 
diagonal elements ^ Vi^Viv We have shown that if 
the rank of pi is £, for any p > £, the operator aju a v docs 
not correlated. Therefore, the potential order operator 
can be constructed as a superposition of nonzero diagonal 
elements of pi 

°i = ^2 W » a U a iV ( 21 ) 

Clearly, the definition still makes the order parameter 
not unique especially if the rank of pi is large. To have 
an explicit form, we can apply various conditions on the 
set of coefficient {w^}. 

The first condition can be the traceless condition, i.e., 

"WV = 0, (22) 

which makes 

tr(Oi) = 0. (23) 

A special case is that the rank of pi is 1, then the traceless 
condition will not leads to any order parameter. With the 
above condition, we need only to calculate the function 
(OiOj) for the correlation function. 

The second condition can be the normalization con- 
dition. That is the maximum value in the set {w^} is 
1. This condition defines a maximum amplitude of the 
correlation function (OiOj). 

Usually, the above two conditions are enough to de- 
cide the form the order parameter. In case that there 
arc still some degrees of freedom in {w^}, one may fix 
the residue degrees of freedom by considering the other 
related physical properties of the ground state. 

Once the order parameter is obtained, the correlation 
function becomes 

(OiO-j) = y^iOilipi^jpj^iifiin'ipjv'lOj) 

/i v 
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^ w^w v V 



/ J 1 
\x'v' 



= ^ WnWvPp'Pv'Pp 'v',,iv (24) 

Clearly, if -P^v,^ is a unity matrix, p i{ jj = pi ® pj, 
(OiOj) — 0, so there is no correlation. Otherwise, 



(OiOA-iOMOJtO, 



(25) 



then Oi plays a role of order parameter. 

Off-diagonal long-range order: the off-diagonal long- 
range order exists if the combined block p iU j is not di- 
agonal in the basis \ipiptpjv) ■ For the block i, since we 

already have (fflLoji/) = ( a L a v) = f° r I 1 ^ v i the 
nonzero off-diagonal elements in pujj means the existence 
of the off-diagonal long-range order. To find the order 
parameter, the first step is to find all possible indices 

pair (p,v) (p 7^ v) with which ^aj^a^aj^Ojv^ ^ and 
a \v a ip. a \^ a jv/ 7^ 0. Then the order parameter can be 



expressed as 

Oi = V (w^ v a\^a lu + w^aL«v) ■ ( 26 ) 



Here we suppose the order parameter is a Hcrmition op- 
erator. From the condition, we have tr(Oj) = 0. We have 

noticed (^al^a^a^aj^ ^ and (al^a^a^ajj^ ^ 0, so 

the correlation function (OiOj) is nonzero too. 

The mode of the order parameter: Having obtained 
the exact form of the order parameter, the next step is 
to determine the mode of the corresponding correlation 
function. For this purpose, we reexpress (OiOj) as a 
function of r = i — j. Usually, the correlation function 
is an oscillating function. If the wavelength of the cor- 
relation function is A. the mode of the order parameter 
is 



k - A • 



(27) 



Finally, we summarize our scheme in the following: 

i) Via studying the mutual information between two 
blocks for various block size, find the minimum block 
size with which the mutual information does not vanish 
at a long distance. 

ii) Calcuate the spectra {p^} of the reduced denisity 
matrix pi of the block i. 

iii) Define the potential diagonal order parameter as 



(28) 



where £ is the rank of the reduced density matrix, and 
a ii |0) defines the eigenvector of the matrix with the 
eigenvalue p^. Apply the traceless condition, normal- 
ization condition, and other suitable restrictions, to de- 
termine the set {w^}. 



iv) If the reduced density matrix of the combined block 
p iU j is not diagonal in the basis of pi ® pj , define the off- 
diagonal order parameter 



+ w*a]„a. 



(29) 



where the pairs (p, v) can be determined from nonzero 
off-diagonal elements in pujj. Determine the set {w^} 
with suitable conditions. 

v) Calculate the correlation function (OiOj) as a func- 
tion ofi—j to determine the mode of the order parameter. 



A. Application I: Ferromagnetic long-range order 

As a simple application, we take a spin chain with fer- 
romagnetic long-range order as an example. We consider 
such a state, 



l*> = ^(ltt 



tt> + iu • • • u» • 



For this state, we can find the minimum block size with 
which the mutual information does not vanish is 1. The 
reduced density matrix of a single spin is 



_ / 1/2 
Pi ~ I 1/2 



(30) 



in the basis of {|f) , ||)}- For any two spins at site i and 
j, the reduced density matrix p iU j 



PiUj 



/ 1/2 \ 





\ 1/2/ 



(31) 



in the basis of {|ft) , III) , lit) , III)}- Therefore, 

S( Pi ) = l, (32) 

and 

S(p iuj ) = l. (33) 

The mutual information between any two spins becomes 

S(i\j) = S ( Pl ) + S ( Pj ) - S ( PlU] ) (34) 
= 1. (35) 

Therefore, there must exist long-range order in the state. 

For the present state, pi is already diagonal and its 
spectra are {1/2,1/2}. So we can define the order pa- 
rameter as 



Oi = wia^an + w 2 a\ 2 ai 2 . 
Then under the traceless condition 



1 



-Wi 



:W 2 = 0, 



(36) 



(37) 
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we have 

wi = —W2- (38) 

Let w\ = 1, we find a traceless operator for the local site 

Oi = a^au - a\ 2 ai2, (39) 

which actually is the Jordan-Schwinger representation of 
the z-component of Pauli matrix of . 

Having found the order parameter, the next step is 
to determine the mode of the corresponding correlation 
function. From the two-site reduced density matrix pmj , 
we have 



(a] 2 a\ 2 a a a, ]2 ) = -, 



(40) 
(41) 



which is independent of the distance between i and j. 
Therefore, the correlation function is 



kM) - <of>W) = i. 



(42) 



So in this case, the order parameter is a z with mode 0. 

For a comparision, we consider another state with anti- 
ferromagnetic long-range order, 

i*) = ^(m---n> + iit---it)). 

In a similar way, we can find that the order parameter is 
a z . Neverthelesse, we notice that for even \i — j\ 



1/2 



I 

1/2 



(43) 



while for odd \i — j\ 



PiUj 





1/2 

1/2 





So the correlation function becomes 



(oftr?) — 1, for even |i — j\ 



T, for odd \i — j\ 



(44) 



(45) 
(46) 



In this case, the order parameter is still a z , however, the 
mode becomes tt. 



where 

[i,j} = ^(\t)i\ih-\i)i\tW. 

The inner product (ip2\ipi) decays expontially with the 
system size, so the state we consider can be written as 



|*> = ^=(|V>1> + |V>2». 



(48) 



The state is the ground state of the Majumdar-Ghosh 
mo del [Tol| . For a single spin, the reduced density matrix 
is 



Pi = 



1/2 
1/2 



(49) 



For two spins separated by a long distance 



PiUj 



/ 1/4 
1/4 
1/4 

\ 1/4, 



(50) 



Therefore pujj — Pi <£> Pj, the S(i\j) — 0. There is no 
long-range spin-spin correlation. The results mean that 
the block including only one spin is too small to find the 
long-range correlation. So we let the block include two 
neighboring spins. The reduced density matrix of two 
spins at site i and i + 1 is 



Pi = 



1/8 

3/8 -1/4 

-1/4 3/8 

1/8 



(51) 



The matrix can be diagonalized as 

Pi 



5 1 

g \<Pifl) (fifi\ + g |Vt,l) 



where 



+ g \<Pi,2) (<PiM + g |^,3) (<Pi,a\ 



\<Pi,a) = -j= {WiU+i) - 



\<Pi,i) = Ititi+i). 
IVi,a) = ^(It4i+i) 
\<Pi,3) = \UU+i) ■ 



(52) 

(53) 
(54) 
(55) 
(56) 



B. Application II: Dimer order 

Now let us consider a uniformly weighted superposition 
of the two nearest-neighbor valence bond state, 

\1>i) = [l,2][3,4]---[i-l,L] 

|^) = [i,l][2,3]-..[L-2,L-l] (47) 



are spin singlet state and spin triplet states respectively. 
Then the entropy 

5 

S (pi) = 3 - - log 2 5. 

For any two pairs of spin, if the distance \i — j\ between 
them is even, the density matrix is of dim 16 x 16. In the 
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basis of spin singlet and spin triplet states, the reduced 
density matrix p iL) j is diagonal 

Piuj = ^ \<Pi,o<Pj,o) (<Pi,o<Pj,o\ + ^I- (57) 
The entropy is 

17 

S(piuj) = 5 - — log 2 17. 

The mutual information between the two blocks can be 
evaluated as 

S(i\j) ~ 0.269. 

The nonzero value of the mutual information between 
two spin pairs means that we can find the potential order 
parameter in their spectra. For this purpose, we define 
it as 

Oi = Wq \ip ij0 ) (<Pifi\ +Wx \ip it i) (y>i,l| 

+w 2 \<Pi,i) (<Pi,2\ + W3 (<Pi,3\ ■ (58) 

Due to the symmetry among the three states 
{\ifi) , |</52) j 1^3) }, we let w\ = u'2 = W3. Then the trace- 
less condition becomes 



Hamiltonian reads 



JV 



5 3 

-wo + —wi = 0. 



(59) 



Upon the normalization condition, we let wq = — 1, we 
have 

5 

Oi = -\(Pi,o){fi,o\ + -^\tpi,i){fi,i\ 

5 5 

+ 3 1^,2) (<Pi,2\ + 3 1^,3) (<fi t a\ (60) 

In terms of spin operators, the order parameter can be 
expressed as 



Oi = 1 + -Oj • <7j +1 . 

The correlation function then can be calculated as 



(O l O J ) = -, for even 
(OiOj) = J, for odd 



(61) 

(62) 
(63) 



The mode of the order parameter is therefore it. 

Traditionally, people use the operator cr, • (Ti+i [ll| to 
describe the spin dimer order. Except for a difference in 
a constant, the order parameter Eq. (|61l) obtained via 
our scheme is exactly the same as the traditional order 
parameter. 



C. Application III: Diagonal versus off-diagonal 
long-range correlation 

To see the off-diagonal long-range correlation, we take 
the ground state of the one-dimensional antiferromag- 
netic Heisenberg model as an example. The model 



H = E ( 5 i S *+i + S * S hi + S i S Ui) > (64) 
j'=i 

where Sj, and Sj are spin-1/2 operators at site j. 
The model has a global SU(2) symmetry, and can be ex- 
actly solved by the Bethe-ansatz method [12J. From the 
Bethe-ansatz solution, the energy spectra of the system, 
its ground state properties, and the thermodynamics of 
the system have been studied explicitly. Moreover, for 
the Heisenberg model, the function (SjSj+ r ) now can be 
calculated by using the generation function method [l3[ . 
However, it is not our motivation to readdress all these 
interesting issue. Instead, to illustrate our scheme, let us 
forget all of the known results at the beginning. 

Now suppose we do not know any information of long- 
range correlation existing in the ground state of the sys- 
tem, we can only calculate the ground state via numer- 
ical methods such as exact diagonalization and density- 
matrix renormalization group[lJ|. With this state, we 
find that the reduced density matrix of a single site is 



_ 1/2 
Pi ~ { 1/2 



(65) 



so the entropy S(pi) = 1. For two separated spins, the 
reduced density matrix has the form 



PiUj 



(u 

w z 

z w 

Vo u 



(66) 



with 2u + 2w = 1, w 5^ it, and z / 0. For instance, 
if |i — i\ = 5, we have u ~ 0.219, w ~ 0.281, and z ~ 
0.062[15j, then S(piuj) — 1.969. The mutual information 
becomes 0.031. Differ from the usual long-range order 
for which the mutual information tends to a constant 
as the distance between two blocks tends to infinite, we 
find the mutual information here decays algebraically and 
becomes zero finally. However, such a correlation still has 
a divergent correlation length, so it still differs also from 
the disordered state for which the correlation length is 
finite. To distinguish the correlation from disorder and 
long-range order, we call such a correlation as long-range 
correlation instead of long-range order. 

Since w 7^ u, according to our scheme, we can find 
that, as discussed in the subsection III(A), the operator 
with diagonal long-range correlation is just of = a\ x an — 
a\ 2 a l2 . 

However, from the expression (|66p . we notice the off- 
diagonal element z is nonzero too. So we predict that 
there exist some kinds of off-diagonal long-range correla- 
tion. The potential operator can be defined as 



wa\ 1 a l2 + w*a\ 2 an 

xO?-yOV. (67) 
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Due to the exchange symmetry between i and j, one can 
prove that (Of Of) = 0. From the expression we 
have 



(^a^a^a^a^ = 0, 



a\iai2a] 2 a ji ) = z -, 



(a^aaa^a^ = 0. 



(68) 
(69) 
(70) 
(71) 



These equalities hold true if (Of Of) = (OfOj). There- 
fore, both Of and Of (or their linear combination) can 
be regarded as the operator which has off-diagonal long- 
range correlation. From Eq. (|57)l . we can see that Of 
and Of are the Jordan-Schwinger representation of the 
x and y components of Pauli matrices erf and erf , respec- 
tively. A careful scrutiny may find that u = w — z, this 
equality which can be verified numerically means that 
(of of) = (of of. 

Having obtained the operators for diagonal and off- 
diagonal long-range correlation, we can find the mode 
of these operators by calculating the corresponding cor- 
relation functions. Take the (ofof +r ) as the example, 
we have (ofof +1 ) = -0.59084, (ofof +2 ) = 0.242716, 
(of of +3 ) = -0.200992, (ofof +4 ) = 0.138610 [H-, there- 
fore the mode of the correlation function (of of +r ) is tt. 
In the same way, the mode of the operators of and of 
are 7r too. 

Up to now, we have obtained the desired operators 
with antiferromagnetic long-range correlations. This 



kind of correlation is consistent with the definition of 
the Hamiltonian. The coefficient of each interaction 
term in the Hamiltonian is positive. The three terms 
SJSj +1 , SjSj +1 , SjSj +1 favor to form antiferromagnetic 
order in the ground state. From our analysis, we have 
(of of) = (of of = (of of). The result can be explained 
from the global SU(2) symmetry of the Hamiltonian. 



IV. SUMMARY 

In summary, we have shown that the non- vanishing be- 
havior of the mutual information at a long-range distance 
means the existence of long-range order. While the mu- 
tual information is operator-independent, we could still 
find the potential diagonal or off-diagonal order parame- 
ter from its spectra. A possible scheme to construct the 
order parameter was provided. Our scheme, as shown by 
also three simple examples, can not only find the diagonal 
order parameter, but also the off-diagonal order parame- 
ter. To find a correct order parameter is very important 
in various studies in the condensed matter physics, our 
scheme is therefore instructive for many condensed mat- 
ter theorists to explore new physics in unknown quantum 
many-body systems. 
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